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Abstract 

We have measured branching fractions of hadronic r decays involving an n meson using 485 
fb -1 of data collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 
We obtain the following branching fractions: £>(r~ — ► K~nv T ) = (1.62 ± 0.05 ± 0.09) x 10~ 4 , 
B(t~ -» X-ttVt) = (4.7 ± 1.1 ± 0.4) x 10~ 5 , B(r~ -» ir~TT°rjiy T ) = (1.39 ± 0.03 ± 0.07) x 10~ 3 , 
and B(t~ — ► A'*~r]z/ T ) = (1.13 =t 0.19 ± 0.07) x 10~ 4 improving the accuracy compared to the best 
previous measurements by factors of six, eight, four and four, respectively. 

PACS numbers: 
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INTRODUCTION 



Hadronic decays of r lepton provide a useful tool for studying QCD phenomena at low 
energy Various decay modes including rj meson(s) are interesting for testing the Wess- 
Ziumino- Witten (WZW) anomaly [J, if, chiral theory @, 3, and relations to e + e cross 
sections following from the conservation of the vector current (CVC) 

We measure the branching fractions of t~ — > K~r]v T (unless specified otherwise, charge 
conjugate decays are implied throughout the paper), K~-k q t\v t , and ix~~ix Q riv T decays, and 
that of r _ — > K ~ (892)r]u T ; the latter is evaluated from the corresponding K~TT r]u T mea- 
surement. Studies of these modes have been reported by CLEO jy, Q, [1] and ALEPH j^j, 
however, most of the results are based on rather low statistics, which do not allow one to 
discriminate between different theoretical predictions. We use a data sample with an inte- 
grated luminosity of 485 fb _1 , corresponding to production of 430 million r-pairs collected 
with the Belle detector at the KEKB asymmetric-energy e + e~ collider flo| . 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect K\ mesons and identify muons 
(KLM). The detector is described in detail elsewhere Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector were used for 
the first sample of 144 fb _1 , while a 1.5 cm radius beampipe, a 4-layer silicon detector and 



a small-cell inner drift chamber were used to record the remaining 341 fb |l2j. 

In this analysis, we use a data sample 100 times larger than any of the previous mea- 
surements. In addition, peaking backgrounds are estimated precisely to decrease systematic 
uncertainties. 

A study of resonance formation in the hadronic final states is in progress and will be 
reported later. 



EVENT SELECTION 

Candidate e + e~ — > r + r~ events are selected with the following common properties: 



r si ^Xrju T and r t + -»• (e//i)+^i/ T , 



where X denotes K~, K~ir°, or 7r~7r°. Candidate rj mesons are reconstructed through two 
decay modes: 77 with a branching fraction of 39.39 ± 0.24% or 7r + 7r~7r° with a branching 
fraction of 22.68 ± 0.35% with 7r° — > 77. In order to remove qq contamination, the tag-side 
t is required to decay leptonically i.e. r~ — > (rv T v\ (£ = el u) (with a branching fraction of 



35.20 ± 0.07%). All branching fractions are taken from Ref. [13 
t~ — * K~r)v T selection 

A candidate event is required to contain either two charged tracks with zero net charge 
and at least two 7's in the 77 — > 77 case, or four charged tracks with zero net charge and at 
least two 7's in the 77 — > 7r + 7r~7r° (tt — > 77) case. A charged track should have transverse 
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FIG. 1: M 77 distributions for K~r]v T selection in (a) rj — ► 77 and (b) rj — > tt~it + it decays. Data 
are fit with a Crystal Ball function plus a second-order polynomial for the background (BG). The 
best fit result is indicated by the solid curve with the BG shown by the dashed curve. 

momentum, p t > 0.1 GeV/c, and —0.866 < cos 6 < 0.956 (8 is the polar angle within the 
detector aperture) where p t and 9 are measured relative to the direction opposite to that of 
the incident e + beam in the laboratory frame, and the 7 should have energy E 1 > 0.05 GeV 
within the same polar angle fiducial region as above. 

To select a sample of r-pairs, the thrust vector and total energy in the center-of-mass 
system are required to satisfy |Vthrust| > 0.8 and 3.0 < .E^tai < 10-0 GeV, respectively. Each 
event is then subdivided perpendicularly to the thrust axis into two hemispheres, the signal 
and the tag. To reduce qq contamination the effective mass of the particles on the tag side 
should satisfy the requirements, M tag < m T (1.78 GeV/c 2 ), while on the signal side the 
mass, M sig , should satisfy 0.70 GeV/c < M sig < m T . 

To reconstruct an 77 from 27's, two 7's with E-y > 0.2 GeV are required in the barrel 
region (-0.624 < cos 9 < 0.833) with at most one 7 with 0.05 < £ 7 < 0.2 GeV in the 
endcap region to allow for initial state radiation. In order to reduce incorrect combinations 
with a 7 from a 7r° decay (denoted hereafter as 7^0), the ^-candidate 7 (7^) should not form 
a 7T° mass with any other 7 (tt° veto selection). In this analysis the 7r° mass window is 
denned as 0.105 GeV/c 2 < M 77 < 0.165 GeV/c 2 (a ±3<r range in the detector resolution). 
Particle identification (PID) uses a likelihood ratio, V x , for a charged particle of species x 
(x = fi, e, K, or ir). V x is defined as V x = L x /H y L y (the sum runs over the relevant particle 
species), where L x is a likelihood based on the energy deposit and shower shape in the ECL, 
the momentum and dE/dX measured by the CDC, the particle range in the KLM, the light 
yield in the ACC, and particle's time-of-flight from the TOF counter [j~4| . For the track 
on the tag side, V e > 0.8 for an e candidate and > 0.8 for a fi candidate with p > 0.7 
GeV/c. On the other hand, a kaon is identified as a track satisfying not only Vk > 0.9, but 
also V e < 0.2 in order to suppress two-photon events, with p > 0.3 GeV/c. 

We also include correlations between tracks and 7's. The polar angle of the missing mo- 
mentum that must be attributed to neutrinos, should satisfy —0.866 < cos 0(P miss ) < 0.956. 
The opening angle between the K~ and 77 satisfies the requirement, cos#(P^ M ; P^ M ) > 0.8. 
The opening angle and energy of two 7^'s should satisfy the following condition: 0.5 < 
cos^(P 7 c 1 M ; P 7 C 2 M ) < 0.96 and (£° M - E™)/(E™ + E™) < 0.8, respectively. 

The 77 mass distribution obtained after these requirements is shown in Fig. 1 (a). 
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FIG. 2: Distribution of likelihood ratios for (a) K~ir r]v T and (b) 7r~-K°r]v T candidates. The points 
with error bars are the data. The yellow, green hatched, and blue histograms indicate the signal, rr 
background, and qq background MC distributions, respectively. The MC histograms are normalized 
to the integrated luminosity of the data. For demonstration purposes, the branching fractions of 
signal MC decays are set to 10 -3 and 10 -2 for (a) and (b), respectively. The dashed vertical line 
show the likelihood ratio requirement. 

In the case of r\ — > 7r + 7r~7r° reconstruction, candidate events should have two additional 
charged tracks compared to the i] — ► 77 sample, but two of the 7's have to form a 7r° instead 
of an j]. While this mode provides one more constraint compared to those in the previous case 
improving the background rejection, the higher multiplicity reduces the detection efficiency. 

The selection criteria different from those in the previous case are indicated below. tt + 
and 7t~ candidates are required to have 0.3 GeV/c < p < 2.0 GeV/c and to be inconsistent 
with the e hypothesis (V e < 0.2) to reject the two-photon background. The photons used to 
form 7r°-candidates are required to have E 1 > 0.1 GeV. In order to remove the contributions 
from higher mass states such as decays of an u meson, the condition M^+^-^o < 0.7 GeV/c 2 
is imposed. The condition on the polar angle of the missing momentum is the same as in 
the previous case while the 7r° momentum should be P^ M > 0.5 GeV/c. 

The resulting 7r + 7r _ 7r° mass distribution is shown in Fig.l (b). 

Selection of r~ — >• K~Tr°r)i/ T and 7r~7r°?7z/ T 

For these r decay modes, we use only 77 —>■ 77 to reconstruct the 77, because of the small 
detection efficiency for 77 — > 7r + 7r _ 7r°. Correspondingly, a candidate event should contain 
two charged tracks and at least four 7's. 

The selection criteria different from those for K~r\v r are listed below. The total momen- 
tum on the signal side is required to satisfy S-P s ^ g M > 2.5 GeV/c; two additional 7's are 
required to lie in the barrel region on the signal side that form the 7r° mass; a condition on 
cosine of the opening angle between the missing momentum and the direction of the thrust 
axis pointing to the signal side is imposed: cos 9(P™ S ; -P t hrust) < — 0.6. For the ^~t^t]v t 
mode, 7T~ candidates should have Vk < 0.2 and V e < 0.2. 

To further suppress backgrounds, we apply a likelihood selection using seven vari- 
ables, such as iKhmstl, P™, ^miss (missing mass squared), P^ M , E° M , E-P s ? g M ; an d 
cos#(P^.; P^ M )- The rr background MC is used for the background likelihood while the 
distributions from K~ir rjv and 7T~7r°r]v MC are used for the signal likelihood in each case. 
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FIG. 3: M 77 distributions for (a) K~ir°r]v T and (b) ir~7r°r]v T candidates. Data are fit with a 
Crystal Ball function plus a second-order polynomial for the background (BG). The result of the 
best fit is indicated by the solid curve with the BG shown by the dashed curve. 



The resulting likelihood ratios denned for K~7r°r]u T and 7r~TT°r]i> T are shown in Figs. 2 (a) 
and (b), respectively. About half of the background is removed while 93% and 90% of the 
signal is retained for each of the respective modes. 

The obtained M 77 distributions are shown in Figs. 3 (a) and (b) for the K~iPr\v T and 
7i^7i°r]i> T modes, respectively. 



BRANCHING FRACTIONS FOR K~rju T , K~TT°r]u T AND vr^vr ^ DECAYS 



To determine the signal yields we determine the number of r/'s and then subtract cross- 
feeds. In order to extract the number of r£s from the resulting M 77 distributions, Figs. 1 



and 3, we use the Crystal Ball function [15j to represent the r] signal and a second-order 



polynomial for the background distribution. The result of the fits is indicated by the solid 
curve in the corresponding figures. The best fits give an 77 mass, m v , of 0.545±0.001 GeV/c 2 , 
in good agreement with the PDG value of 0.54751 ± 0.00018 GeV/c 2 , and a resolution 
o m = 0.012 ± 0.002 GeV/c 2 for the three different decay modes in the 77 — > 77 case and 
mj= 0.5474 ± 0.0007 GeV/c 2 with a mji = 0.0075 ± 0.0004 GeV/c 2 for the K~r]v T mode in 
the 77 — > 7r + 7r~7T° case. 

The 77 yields obtained from the fits are Nk-t) Vt — 1387 ± 43, N^-^o^ = 270 ± 33, and 
N n - n Q„ VT = 5959 ± 105 events for the K~r]u T , K~it°riv T , and ir~it°riv T modes, respectively, 



in the 77 — > 77 case. The yield for the 77 — > 7r + 7r _ 7r° case is Nx~ V u T -,ri-*3n — 241 ± 21 events. 
The rj signal region is defined as 0.48 GeV/c 2 < M 77 < 0.58 GeV/c 2 in the 77 — > 77 
case while it is 0.52 GeV/c 2 < M^+^-^o < 0.58 GeV/c 2 in the 77 — >• 7r + 7r~7r° case. 
The detection efficiencies are estimated with MC simulation in the same manner as the 
detection of 77 yields. The corresponding efficiencies are e = 0.94 %, 0.35 %, 0.47 %, and 
0.16 %, respectively, and include the intermediate branching fractions such as B(rj — > 77), 
B(r] — > 7r + 7T~7r°), and B{t~ — > l~v T i>\). 

These event yields include backgrounds classified into three categories. One is due to 
cross-feed effects between the three signal modes, which is taken into account by solving the 
following system of equations: 

N k - Wt = 2N TT (B(K-T)u T ) ■ e{ + B(K-Tr°r)v T ) ■ e\ + B{tx~tx\v t ) ■ el), (1) 
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TABLE I: Raw r\ yields for the four selected r decay modes. N„ is the total number of r\ events 
detected, which include cross-feed from two other modes. The contributions of qq and 'other', 
mostly 7r _ 7r 7r 77^ r and ir~ir + Tr~r]i' T , are also included. 



Mode 


N v 




Kr)v T 




TT-K 7]U T 


Other 


qq 


K~j]v T (77 -> 77) 


1387 ± 43 






15.1 ±3.8 


18.0 ± 1.0 


1.1 ±0.2 


30.6 ± 15.6 


K~ir G rju T 


270 ± 33 




16.0 ±0.9 




85.3 ±4.6 


1.2 ±0.4 


27.0 ±8.5 


TT^TT°r]lJ T 


5959 ± 105 




2.4 ±0.1 


9.4 ±2.4 




71.6 ±20 


212 ±29 


K~r]v T {rj — ► 7r + 7r _ 7r°) 


241 ± 21 






3.3 ±0.8 


5.8 ± 1.3 


< 1.18 


9.1 ±2.2 



Nk-^ = 2N TT (B(K-r]u T ) ■ e\ + B{K~^ Wt ) ■ e\ + S^Vr) • el), (2) 

N^o VUt = 2N tt (B{K- Wt ) ■ e\ + B{K~ix\v T ) ■ e\ + B(tT^Vt) • ejj), (3) 

where B(K~r)v T ), B(K~7r°r]i> T ), and B(7T~7T r]u T ) are the branching fractions of the respective 
modes and N TT is the total number of r pairs produced. Here is the detection efficiency 
in each case and j (i) indicates the decay sample (selection criteria). 

The second type of background is from decay modes of the r-lepton itself, such as 
7r _ 7r 7r°?7z/ T and 7T~7T + 7T~r]i> T . These backgrounds are estimated by using a rr MC simu- 
lation with branching fractions taken from [l6[. The estimated background is included in 
the 'other' category in Table [H The contribution of this type of background is negligible 
for the K~7]v T and K~ir°riv T modes, and is smaller than the statistical uncertainty in the r\ 
yield obtained above for the ir~it°riv T mode. The contamination from n~riv T decay should 
also be considered. This decay is strongly suppressed since it violates G-parity and proceeds 
via a second-class current. Its branching fraction is predicted to be 10 -5 . Therefore, its 
contribution to each mode is negligible. 

The last background category is e + e~ — > qq, which is estimated from MC simulation. The 
MC was tuned beforehand and validated with a qq enriched sample, which was produced with 
some variations of the event selection criteria. The M S i g cut is removed and the condition 
M tag > m T is implemented on the tag side. In addition, the PID requirement for the charged 
track on the tag side was reversed (i.e. V e < 0.8 and < 0.8). The rj yield of MC was then 
tuned to be consistent with that of data. The qq contributions are 2-4 % of the rj yields for 
the K~7]v T and it~ir°riv T modes and 10% for the K~n°riv T . They also are summarized in 
Tabled 

After subtracting the 'other' and qq contributions from the individual r\ yields, we solve 
the system of equations to obtain the branching fractions for three decay modes. They are 
(1.62 ± 0.05) x 1(T 4 , (4.7 ± 1.1) x 10" 5 , and (1.39 ± 0.03) x 1(T 3 for K~r]u T , K~n Q r)v T , 
and 7r~7r°?7^ T , respectively. The cross-feed yields obtained are also listed in Table [B The 
number of cross-feed events for K~j]v T in the 77 — > 7r + 7r~7r° case is evaluated using the above 
branching fractions, obtained in the 77 — > 77 case. The branching fraction for K~7]v T in this 
case is (1.65 ±0.16) x 10~ 4 . 

The systematic uncertainties are estimated as follows: the estimation of peaking back- 
grounds provides sizable uncertainties only in case of the 7i~7r°riis T (3.3 %) and qq (6.0 %) 
contributions to the K~^i]v T mode. As for the qq estimation, due to the finite statistics of 
qq enriched sample, the uncertainties of 26 %, 19 %, and 9.6 % for K~t]u t , K~-k^t\v t , and 
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TABLE II: Summary of systematic uncertainties in each mode (%) 



Mode 


K T)V T 


K -K Q rju T 


7T TT°r]U T 


K~T]V T 


K*-r)v T 


r\ detection 


77 -» 77 


T] -> 77 


T] -> 77 


77 — ► 7r + 7r _ 7r 


77 -> 77 


Estimation of K~rju T 


— 


0.6 


1.8 x 10~ 3 


— 


— 


Estimation of K~ifir\v T 


0.3 


— 


4.2 x 10~ 2 


0.4 


— 


Estimation of TT~TT rju T 


7.5 x 10~ 2 


3.3 


— 


0.1 


— 


Estimation of 7r _ 7r 7r°Ty^ r 


— 


— 


0.4 


— 


— 


Estimation of qq 


1.5 


6.0 


0.5 


1.5 


2.4 


Particle ID (K/ir) 


3.3 


2.2 


1.0 


2.8 


2.2 


Particle ID (Lepton) 


2.3 


2.8 


2.6 


2.6 


2.6 


Track finding 


1.3 


1.3 


1.3 


3.3 


1.3 


Luminosity measurement 


1.6 


1.6 


1.6 


1.6 


1.6 


7T° detection 




2.0 


2.0 


2.0 


2.0 


7T° veto 


2.8 


2.8 


2.8 




2.8 


Signal MC 


0.5 


1.7 


0.7 


1.3 


1.7 


B(rj — » 7T + 7r~7r ) 








1.6 




Total 


5.6 


8.9 


5.0 


6.3 


6.0 



r K~ir 'qv T decays arise from tuning, respectively. The errors in the qq background estimations 
in Tabled] come from this uncertainty and the statistical uncertainty in the qq MC, and are 
treated as systematic uncertainties. The uncertainties in the peaking backgrounds in all 
other cases are rather small. Uncertainties in the PID efficiency and fake rate are evaluated 
to be 2-3 % for kaon ID, 1 % for 7r ID and around 2.5 % for the lepton ID; these values 
are obtained by averaging the estimated uncertainties depending on momentum and polar 
angle of each charged track. For the ir° veto selection, the efficiency was compared between 
data and MC with a sample in which the PID of the charged track on the tag side was 
reversed from the usual TT~ir Q riv T selection. The efficiencies were consistent. Therefore, 2.8 
% of the statistical uncertainty from this comparison was counted as systematic uncertainty. 
Other systematic uncertainties are summarized in Table [Til The total systematic uncertain- 
ties are 5.6 %, 8.9 %, 5.0 %, and 6.3 % for the K~r\v T (77 — > 77), K~7r rju T , 7i~7i°r]b> T , and 
K~7]v T (77 — > 7r + 7r~7T°) modes, respectively. 

Taking the systematic uncertainties into account we obtain the following branching frac- 
tions: 

(1.62 ± 0.05 ± 0.09) x 10" 4 for 77 -> 77, 
(1.65 ± 0.16 ± 0.10) x 10~ 4 for 77 -> tt+tt^tt , 
(4.7 ± 1.1 ±0.4) x 10" 5 , 
(1.39 ±0.03 ±0.07) x 10" 3 . 

By combining two measurements for r~ — > K~r]v T decay, we obtain: 

B{t~ -> K-r]u T ) = (1.62 ± 0.05 ± 0.09) x 10~ 4 . 



B(t~ -> K~r]v T ) = 
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M^o , GeV/c 2 

FIG. 4: The K~tt° invariant mass distribution for the K~irr]v T events. Data are fitted with 
a Breit-Wigner (BW) function in (a), and the best fit is indicated by the solid curve while the 
continuum component evaluated from sideband regions is shown by the dashed curve. The fit 
gives N k *- VUt = 119 ± 19 events with x 2 /(d.o.f. = 27) = 1.15 (the probability to obtain this 
result is 0.265). (b) and (c) show the results of similar fit with a phase space distribution only 
and a BW plus a phase space distribution. In these cases the results are N^*--^ = 102 ± 21 
events with x 2 /(d.o.f. = 27) = 2.09 (the probability is 0.0008) for the phase-space distribution, or 
N K *-r,u T = 105 ± 20 and N K -^ VUr = 33 ± 30 events with x 2 /(d.o.f. = 26) = 1.13 (the probability 
is 0.294). 

BRANCHING FRACTION FOR K*~ (892)t/i/ t 

From the K~n®riv T samples within the r\ mass region, 
0.50 GeV/c 2 < M 77 < 0.58 GeV/c 2 , we extract a branching fraction for t~ — > 
K*-(892)r]u T decay through the if*" (892) -> K~n° mode (B(K*~ (892) -> K~n°) = 1/3). 

The distribution of K~ir° invariant mass, M^-^o, for the selected samples is shown in 
Fig. |H The r~ —* 7i~ir riu T background, which cannot be neglected as shown in Table 
[H is estimated by using MC simulation with the branching fraction measured in this pa- 
per. The other backgrounds are estimated from two sidebands of the M 77 distribution: 
0.43 GeV/c 2 < M 77 < 0.48 GeV/c 2 and 0.60 GeV/c 2 < M 77 < 0.65 GeV/c 2 . The 
background distribution is indicated by the dashed curve. The peculiar shape of the ex- 
pected background is due to two components contributing to it: that of TT~7r r]u T at high 
mass and the one from all other r decays at low mass. A clear excess above the background 
is seen around 0.9 GeV/c 2 , suggesting if* - (892) resonance formation. Three types of signal 
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TABLE III: Comparison of our measurement with previous results 



H It 1 

Mode 


y~y ■ , l ■ i 

£> m this analysis 


Previous o 


Kelerence 


— > K~T]V T 


(1.62 ±0.05 ±0.09) x 10~ 4 


(2.6 ±0.5 ±0.5) x 10" 4 
(2.9 ± 1.3 ±0.7) x 10" 4 


CLEO [7] 
ALEPH [9J 


T~ — ► K^^'qi'r 


(4.7 ±1.1 ±0.4) x 10~ 5 


(17.7 ±5.6 ±7.1) x 10~ 5 


CLEO [8] 


T~ — > 7r~7r°7/!/ T 


(1.39 ±0.03 ±0.07) x 10~ 3 


(1.7 ±0.2 ±0.2) x 10~ 3 
(1.8 ±0.4 ±0.2) x 10~ 3 


CLEO [6] 
ALEPH [9] 


t~ — ► K*~r]u T 


(1.13 ±0.19 ±0.07) x 10~ 4 


(2.90 ±0.80 ±0.42) x 10~ 4 


CLEO [8] 



functions have been tested. A Breit-Wigner (BW) function whose mass and width are set 
to those of the A'* - (892) is fitted to data in Fig. @] (a). A fit with the hypothesis that the 
excess signal events are due to a non-resonant, but V — A pure phase space process is shown 
in Fig. H] (b). Figure 0] (c) shows the result of a fit with a BW plus a phase space function 
as a signal. In each case, the signal function is smeared to take into account the detector 
resolution. The BW function describes the data well in Fig. H] (a) and the resulting number 
of if*~(892) events is N k *- Wt = 119 ± 19. In contrast, the fit shown in Fig. |4] (b) is in 
complete disagreement with the data. Although inclusion of a small phase space component 
in addition to the dominant K*~ resonance component also represents data well as shown 
in Fig. H] (c), the resulting magnitude of the phase space component is consistent with zero 
within errors. Therefore, we conclude that all excess events are produced via the i^*~(892) 
resonance. 

The detection efficiency is evaluated by MC simulation to be e = 0.12%, including the 
branching fractions of i3(A'*"(892) —> K~n°), E{r\ — > 77), and B(t — > iu T v). 

Since the requirement of A'*" (892) formation in this case is a rather strong constraint, 
no significant peaking background contribution from other r decays is found. Therefore, the 
systematic uncertainty in the A'*~(892)?7z/ r mode is small as compared with the K~^r\v T 
analysis. Other sources of systematic uncertainties are the same as those of K~ir°riv T , except 
for the background contamination; a total systematic uncertainty of 6.0 % is obtained with 
details summarized in Table [Til Finally, we obtain the following branching fraction: 

B{j- K*~T)v T ) = (1.13 ± 0.19 ± 0.07) x 10" 4 . (4) 



RESULTS AND DISCUSSION 



We have obtained branching fractions for four different decay modes based on a high- 
statistics data sample of 430 million r-pairs collected with the Belle detector: 

B(r~ -> K~r]u T ) = (1.62 ± 0.05 ± 0.09) x 10" 4 , 

B(t~ -> ATttVt) = (4.7 ± 1.1 ± 0.4) x 10" 5 , 

B(t~ -»■ TTTrVr) = (1.39 ± 0.03 ± 0.07) x 10" 3 , 

B(r~ -> K*"r]u T ) = (1.13 ± 0.19 ± 0.07) x 10~ 4 , 

where the first error is statistical and the second is systematic. 
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FIG. 5: The invariant mass distributions of various combinations of final state particles for t~ — > 
7r~7r°r]i' T decay. The points with error bars are the data. The yellow, green hatched, and blue 
histograms indicate the signal, rr background, and qq background MC distributions, respectively. 
The dominant backgrounds come from rr events. The qq background is strongly suppressed and 
negligible in our sample. 



Compared to previous experiments, we have improved not only the statistical uncertain- 
ties, but also evaluated reliably the background contamination. In Table HTT1 our results are 
compared to those previously obtained by the CLEO 0, 0, 0] and ALEPH j9| collabora- 
tions. Our measurement improves the uncertainties in the branching fractions by a factor 
of six (K~riv T \ eight [K~^riv T \ four (7r~7r°?7i/ T ), and four (K*~r]u T ) compared to the most 
precise determinations from CLEO. In addition, the high statistics of our experiment allows 
much more reliable estimation of various backgrounds including the peaking one. The rela- 
tively poor statistics of previous measurements imposed some limitations on BG estimations 
(l7| . It is also noteworthy that in all cases the central value of our measurement is lower 
than that of the other measurements. This may be due to the underestimation of back- 
grounds in the previous experiments. The improved accuracy in the branching fractions of 
the decay modes reported here is important for future searches for the second-class current 

— > ix~r\v T decay. 

Our branching fraction for r _ — > -K^-K°r]v T decay is consistent with the prediction based 
on CVC and experimentally measured e + e~ — > 7i + 7i~r] cross sections jH|. In addition, the 
Monte Carlo code TAUOLA reproduces the observed hadronic mass distributions rather 
well as shown in Fig. [S] while more tuning is needed for r's decay modes involving kaon(s). 
The values of the branching fractions obtained for r _ — > K~i]v T and — > K~TT°rjv T decays 
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differ slightly from the predictions by Li [4(. 

Further studies of final state dynamics and resonance formation in the r~ — > K~ir°riv T 
and n~TT°r]u T decays, which may be important for understanding the WZW anomaly are in 
progress. 
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